Correlated-electron description of the photophysics of thin films 

of 7r-conjugated polymers 

D. Psiachos and S. Mazumdar 
Department of Physics, University of Arizona Tucson, AZ 85721 
(Dated: January 20, 2009) 

Abstract 

We extend Mulliken's theory of ground state charge transfer in a donor-acceptor complex to 
excited state charge transfer between pairs of identical 7r-conjugated oligomers, one of which is in 
the optically excited state and the other in the ground state, leading to the formation of a charge- 
transfer exciton. Within our theory, optical absorptions from the charge-transfer exciton should in- 
clude a low energy intermolecular charge-transfer excitation, as well as distinct intramolecular exci- 
tations from both the neutral delocalized exciton component and the Coulombically bound polaron- 
pair component of the charge-transfer exciton. We report high order configuration-interaction cal- 
culations for pairs of oligomers of poly-paraphenylenevinylene (PPV) that go beyond our previous 
single configuration-interaction calculation and find all five excited state absorptions predicted us- 
ing heuristic arguments based on the Mulliken concept. Our calculated excited state absorption 
spectrum exhibits strong qualitative agreement with the complete wavelength-dependent ultrafast 
photoinduced absorption in films of PPV derivatives, suggesting that a significant fraction of the 
photoinduced absorption here is from the charge-transfer exciton. We make detailed comparisons 
to experiments, and a testable experimental prediction. 

PACS numbers: 42.70.Jk, 71.35.-y, 78.20.Bh, 78.30.Jw 
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I. INTRODUCTION 



It is now generally recognized that interchain interactions play a strong role in the pho- 
tophysics of thin films of 7r-conjugated polymers (PCPs).-^^^^^^^^^ 1 ^^ 1 ^ While 
the primary photoexcitation in dilute solutions of PCPs, which contain widely separated 
single strands, is to the optical exciton, interchain species other than the intrachain opti- 
cal exciton are generated in films. Since the original demonstration of interchain species 
in films of poly(2-methoxy,5-(2'-ethyl-hexyloxy)l,4 paraphenylenevinylene) (MEH-PPV) by 
Rothberg and collaborators,- a variety of complex optical phenomena in many other PCPs 
have been ascribed to interchain interactions. Experiments that indicate the formation of 
interchain species upon photoexcitation include, (i) photoluminescence (PL) at longer wave- 
length relative to solutions, with smaller quantum efficiency,-^^^ 1 ^ (ii) delayed PL lasting 
until milliseconds, and its quenching by an electric field, -^^^ and (iii) additional photoin- 
duced absorptions (PAs) in ultrafast spectroscopy of films over and above those observed in 
solutions.-* 1 ^ 1 ^ 

Three different kinds of interchain species can in principle occur in PCPs. These include, 
(i) neutral exciton delocalized over multiple chains, as in an aggregate,— ii&i2j2£ (ii) com- 
pletely ionic polaron-pair, with positively and negatively charged neighboring chains that 
are Coulombically bound,-* 2 -! 4 -, and (iii) the charge-transfer (CT) exciton, which is a quan- 
tum mechanical superposition of the delocalized neutral exciton and the ionic polaron-pair 
that is obtained in the presence of significant interchain electron hopping. 4 ' 21 All three kinds 
of interchain species have been proposed to explain experiments in real systems. PL at 
wavelengths longer than in solutions with weaker quantum efficiency can be explained using 
the H-aggregate scenario, wherein absorption is to the highest state in the exciton band 
and emission occurs from the lowest state in the band.— ^ On the other hand, PA in films 
was originally thought to be from a polaron-pair state.-! 4 - Very recently, Wang et alM have 
claimed that the detailed features of the PA in films 1 ^* 1 ^ can be understood only within 
the CT exciton picture. According to Wang et al. a significant fraction of both PL and 
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were based on 



PA occur from the CT exciton. The calculations of the PA in Reference 
the single configuration-interaction (SCI) approximation, and were able to reproduce only 
a limited low energy region of the experimental PA spectra.— 1 ^ In the present work we go 
beyond SCI and perform quadruple CI and multiple reference singles and double CI (MRS- 
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DCI) calculations^ 1 ^ of excited state absorptions in interacting pairs of PPV oligomers to 
understand the PA spectra over the complete experimental energy range. Our calculated 
PA spectra, originating from the CT exciton, exhibit strong qualitative agreement with the 
experimental PA spectra in films of PPV derivatives,— in spite of the short lengths of the 
oligomers for which the many-body calculations could be performed. 

It is important for what follows that we clearly state the difference between PA in solutions 
and films. Both solutions and films show absorptions labeled PAi and PA2. The lower energy 
PAi exhibits a peak at ~ 1 eV and the higher energy PA 2 occurs at ~ 1.3—1.4 eV. Two distinct 
additional PAs, Pi at ~ 0.4 eV and P2 close to, but distinct from PA2, are seen in films.— ^ 
Whether or not Pi and P2 arise from free polarons (generated by exciton dissociation) or from 
a bound interchain species has been controversial in the past (see discussions in reference I2I ) . 
Recent experiments have found similar behavior in heterostructures consisting of donor and 
acceptor PCPs,— & or of PCPs and finite molecular acceptors.— Identifying the nature of the 
dominant intermolecular species, and determining their photophysics is clearly important. 

In section II we discuss our theoretical model and the methods we use. In view of 
the highly correlated natures of the two-chain wavefunctions, and the complexity of the 
phenomenon of excited state charge transfer that is only beginning to be understood, we 
present in section III heuristic physical descriptions of all photoexcitations that are possible 
from the CT exciton. This discussion is an extension of Mulliken's theory of ground state 
charge transfer between a donor and an acceptor— to excited state charge transfer between 
identical pairs of molecules, one of which is in the optically excited state. The qualitative 
arguments make the computational results, presented in section IV, understandable. We 
present our conclusions, along with comparisons to experiments, in section V. 



II. THEORETICAL MODEL AND METHODS 



As in Reference 
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we consider pairs of capped PPV oligomers 4 nm apart, stacked co- 
facially. The sophisticated many-body approaches used here require enormous Hamiltonian 
bases,— which necessarily limits us to oligomers that are considerably shorter than those in- 
vestigated in our previous SCI calculations. ?1 Our calculations here are for pairs of oligomers 
of length three and four units each. We refer to the individual oligomers as PPV3 and PPV4, 
respectively. Note that with two such chains these calculations are much larger than the 
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previous QCI and MRSDCI calculations of PAi and PA2 in single chains.— Although the 
actual arrangements of the chromophores deviate from ideal cofacial stacking in the real 
systems, it is believed that this assumption captures the essential physics of polymer thin 
films.— 'i^i Unlike in |2l|, our calculations here are only for pairs of oligomers with equal 
lengths, so that the entire two-chain structure possesses a center of inversion. This is for 
computational simplification only and does not affect our conclusions. As indicated in our 
previous work, lack of inversion symmetry makes the CT exciton below the optical exciton 
weakly optically allowed from the ground state, but the energy locations and strengths of 
the calculated PAs originating from this interchain species are the same whether or not 
inversion symmetry exists.— 

We describe our system within an extended two-chain Pariser-Parr-Pople 7r-electron 
Hamiltonian^i^ H = Hi ntra + Hi nter , where Hi ntra and Hi nter correspond to intra- and 
interchain components, respectively. The two components of H are written as, 



fj,(ij),a n,i 



(1) 



and, 



H inter = ~ t ± ( C l,i,T C ^' + H.C.) 

/j.</j,' ,i,a 



(2) 



fJ,<H',i,j 

In the above c i a creates a 7r-electron of spin a on carbon atom i of oligomer //(= 1, 2), 
n ^,i,a = c uia c vAv * s the num ber of electrons on atom i of oligomer // with spin a and 
n Av = Yla n i*,i,a ^ S the total number of electrons on atom % of oligomer \i. We consider 
standard intrachain nearest neighbor one-electron hopping integrals £y = 2.4 eV for phenyl 
C-C bonds, and 2.2 (2.6) eV for single (double) C-C bonds, respectively— U and are the 
on-site and intrachain intersite Coulomb interactions. Vij are obtained from a modification 
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of the Ohno parametrizatio 



where Rij is the distance between carbon atoms i and j in A, and k is the background 
dielectric constant along the chain direction that takes into account screening effects from 



Vij = . = (3) 
k+ 1 + 0M17R% 



the inner-shell electrons. Within the standard Ohno parametrization, U = 11.26 eV and 
K — 1. We use U = 8 eV and k — 2. The detailed justification of our parameters has been 
given elsewhere.— 

As in our previous work,— interchain hopping is restricted between the closest pairs of 
carbon atoms; we choose t± = 0.1 eV. For VA, we choose the same functional forms as 
in Eq. [3j with the same dielectric screening for the PPV oligomer. As shown before,— 
the results of the PA calculations are relatively insensitive to variations of the transverse 
dielectric constant, as long as reasonably realistic values are chosen. 

In addition to PA calculations for interacting pairs of neutral chains, we have also calcu- 
lated the ground state absorption energies of singly- charged three- and four-unit oligomers, 
for comparison to the PA spectra of neutral systems. These absorption energies can be 
strongly affected by electron-phonon interactions, which, however, have not been included 
in our Hi ntra - This is primarily because calculating PA within a model that incorporates both 
electron-electron and electron-phonon interactions is currently difficult, and the comparison 
would be meaningless if only one of them was calculated with electron-phonon interactions. 
We shall nevertheless refer to the charged oligomer absorptions as polaron absorptions. 

III. CHARGE-TRANSFER EXCITONS AND EXCITED STATE ABSORPTION 

We present here a brief review of our earlier work on the CT exciton in interacting PCP 
chains for the sake of completeness. Following this, we present a diagrammatic description of 
the excited state absorptions that are expected for the two-chain system, within an extension 
of the Mulliken theory of charge transfer.— 

Four distinct basis functions, two neutral and two ionic, are relevant for understanding 
the correlated eigenstates near the optical edge of PCPs.— These are the excited neutral 
and charged eigenstates of H intra for two independent chains. We write the neutral two- 
chain basis functions as product functions of the single-chain 1A S ground state and the 
1B U optical exciton, viz., \excl) = \lB u ) 1 \lA g ) 2 and \exc2) = \lB u ) 2 \lA g )x, respectively, 
where the suffixes are chain indices. The charged polaron configurations, in their low- 
est states, are similarly written as \P + )i\P~)2 and \P + )2\P~)i- For symmetrically placed 
chain pairs, and within the SCI, the total Hamiltonian H has four eigenstates:— the even 
parity optical exciton, \EX) = -k=[\excl) + \exc2)]; the even parity polaron-pair state, 
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\PP)+ = -^[\P + )i\P~)2 + |-P + )2|-P~)i]; an d two distinct charge-transfer excitons, here- 
after \CTX1) and \CTX2), that are superpositions of the odd parity states 4=[|excl) - 
|exc2)] and -^[\P + )i\P~)2 — \P + )2\P~)i\- F° r the interchain dielectric screening constant 
less than or equal to the intrachain dielectric screening, the energy ordering of the states^ 
is \CTXl) < \EX) < \PP)+ < \CTX2). 

Within the Mulliken theory— a weakly coupled donor-acceptor system is expected to 
exhibit a low energy charge-transfer absorption as well as weakly perturbed molecular ab- 
sorptions. Current pump-probe experiments have very high sensitivity with \AT/T\ ~ 
(where T is the transmission of the probe beam without the pump beam, and AT the differ- 
ence in the transmission with and without the pump) 1 ^. Under this condition, excited state 
"molecular absorptions" from \CTX1) will be visible from both the neutral exciton compo- 
nent as well as the ionic polaron-pair component. We show this schematically in Fig. 1, 
where we have included all possible excited state absorptions from \CTX1). Here \ExJ) 
and \PP~) are the superpositions of the odd parity exciton and polaron-pair configurations 
mentioned above; we have shown only one member of each symmetry-adapted pair in both 
cases. Each two-chain configuration is shown as a product state of single-particle configu- 
rations, with electrons occupying bonding and antibonding molecular orbitals (MOs), and 
spin-singlet bonds between the MOs that are singly occupied by electrons. 

The charge-transfer excitation from \CTX1) in Fig. 1 is to the even parity polaron- 
pair state \PP) + . This state is reached by intermolecular transfer of an electron from the 
antibonding (bonding) MO of the molecule in the excited (ground) state of \ExJ) to the 
antibonding (bonding) MO of the molecule in the ground state. The polarization of this 
transition, in contrast to the intramolecular transitions discussed below, is perpendicular to 
the plane of the molecules. 

Two different intramolecular excitations from the \ExJ) component of \CTX1) are indi- 
cated in Fig. 1. The first of these is to a higher energy effective one electron-one hole (le-lh) 
excited state that we have labeled as the \mA g J). The second is an effective two electron- 
two hole (2e-2h) excitation labeled \kA g ^). Corresponding to each of these intramolecular 
excitations of the two-chain CT exciton there exist single-chain excited states that have 
been extensively discussed in the context of nonlinear absorptions in single chains.— ^ The 
single-chain \mA g ) and the \kA g ) are the final states corresponding to induced absorptions 
PAi and PA2, respectively, from the single-chain \1B U ) exciton.— ^ The minus subscripts in 
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FIG. 1: Schematic of the PAs from the charge-transfer exciton \CTX1), shown in the box as 
a superposition of the neutral exciton state \Ex-) and the polaron-pair state \PP-) (see text). 
\PP+) is reached by low energy charge-transfer excitation; \mAg_) and \kAg_) are reached by 
le-lh and 2e-2h intramolecular excitations from the \Ex—) component of \CTX1), giving induced 
absorptions PA' X and PA2, respectively. \PP1*) and \PP2*) are reached by le-lh intramolecular 
excitations from the \PP~) component of \CTX1) (see text). 



\mAg_) and the \kA g _) imply that as with the \Ex-), the actual wavefunctions contain odd 
parity superpositions of two-chain configurations. The \mA g _) and the \kA g ^) are higher 
energy CT excitons, and also contain excited polaron-pair components that are not explic- 
itly shown in the Figure. Indeed, each final state wavefunction in Fig. 1 is a superposition 
of several dominant configurations, only one of which is shown explicitly in the Figure. The 
others are obtained by applications of mirror-plane and charge-conjugation symmetry op- 
erations on the configurations shown. In the case of the \mA g J) the exciton and excited 
polaron-pair contributions to the wavefunction were both demonstrated within the SCI in 
our previous work. 21 

Fig. 1 also indicates intramolecular absorptions from the \PP~) component of \CTX1). 
We have labeled the dominant final configurations here as \PP1*) and \PP2*), which are 
bound high energy excited polaron-pair configurations. For realistic t± and V^~, the energies 
of the excited state transitions from the CT exciton to the polaron-pair states \PP1*) and 
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\PP2*) should be comparable to the ground state absorption energies in charged chains, 
since the former are exclusively from the charged component of \CTX1). 

IV. COMPUTATIONAL RESULTS 

We find from direct computations of excited state spectra that all of the absorptions 
indicated in Fig. 1 should indeed be observed. We report results of QCI and MRSDCI 
calculations on pairs of PPV3 and PPV4 oligomers. In our many-body calculations we 
have retained all innermost one-electron delocalized bonding and antibonding band levels, 
and the localized bonding and antibonding bands that occur below and above these.— The 
dimension of the QCI Hamiltonian matrix for pairs of PPV3 oligomers was 1,83,3276. In 
order to obtain convergence in energies to 0.02 eV, the MRSDCI calculations for PPV4 
were performed with 54 reference configurations. The overall dimension of the Hamiltonian 
matrix for the two-chain system in this case was 1,46,8048. These sizes are at the limit of 
our computational capability. 

In Figs. 2(a) and (b) we have shown our calculated PAs from \CTX1) for PPV3 and 
PPV4 oligomers, respectively. The PA energies are given in units of the energy of the two- 
chain optical exciton \EX) in each case. We have also calculated the optical absorption 
energies from the ground states of the charged single chains (PPV3 + and PPV4 + ). The 
arrows in the figures denote the charged chain absorption energies, relative to the energy of 
\EX). The large differences between the charged-chain absorptions and PAs labeled PP1 
and PP2 are primarily due to the finite size effects associated with the very short system 
sizes being probed here, although additional contribution to this difference arises also from 
interchain Coulomb interactions that influence the energies of the two-chain PAs. 

We have calculated the PAs from the two-chain optical exciton \EX). PAs nearly identical 
in energy to PA' X and PA' 2 are found, but the PAs labeled CT, PP1 and PP2 are missing 
in this case, confirming the analysis of Fig. 1. We have also done detailed wavefunction 
analyses of the final states of all the PAs shown in Fig. 2 to confirm that the dominant 
configurations in each case are the same as that indicated in the schematic Fig. 1. Thus 
the final state of the absorption labeled CT in Fig. 2 is predominantly \PP + ), while the 
final states of the absorptions labeled PP1 and PP2 are predominantly superposiitons of the 
polaron-pair configurations shown in Fig. 1. The results of Fig. 2 go considerably beyond 
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the SCI studies of reference 
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which had found only the PA labeled PA[ and the low 
energy charge-transfer absorption. While it is expected that inclusion of higher order CI 
is essential for detecting the PA labeled PA^L, it is interesting that this is a requirement 
also for the PAs labeled PP1 and PP2, in spite of their predominantly le-lh characters. A 
few additional points need to be emphasized. Our earlier SCI calculations were for pairs 
of oligomers 2-3 times longer,— and from careful finite-size analysis we concluded that the 
features labeled CT and PA[ will continue to be seen for a broad range of realistic oligomer 
lengths. While it is not possible to do a similar finite size analysis for the QCI and MRSDCI 
calculations we report here, based on the previous analysis we are convinced that all PAs in 
Fig. 2 will similarly be seen at realistic chainlengths. The relative energy separations may, 
however, be strongly chain-length dependent. Thus from the nature of the final states \PP+) 
and \PP*) in Fig. 1 it is conceivable that the absorptions labeled CT and PP1 in Fig. 2 
are energetically much closer in systems with realistic lengths, especially if electron-phonon 
interactions play a weak role. Similarly, the relative locations of PP2 and PA' 2 may also 
depend on system size. Both of these may have experimental consequences, as discussed in 
the next section. 



V. DISCUSSIONS AND COMPARISON TO EXPERIMENTS 

Within our theory, photoexcitation in films generates predominantly the optical exci- 
ton, along with small amounts of the CT exciton, which is weakly optically allowed due 
to disorder.— PL is predominantly from the CT exciton, whose neutral component is the 
same as the lowest forbidden exciton state in a H-aggregate. The mechanism of light emis- 
sion is therefore the same as in theories emphasizing aggregate formation.— >i2i2£ We believe 
that PA in films at the earliest times is predominantly from the optical exciton with weak 
contributions from the CT exciton, while at later times the PA is predominantly from the 
CT exciton. Recent careful measurements have indicated that the PA labeled PAi in films 
exhibit biexponential decay, which would agree with the two-state scenario.— Recent ex- 
periments on pentacene films support this viewpoint. 35 Additional support comes from the 
absence of free polarons, as noted in microwave conductivity measurements^ and terahertz 
spectroscopy— in spite of the experimentally observed PA bands Pi and P2 that resemble 
absorption from charged chains. We believe that an alternate possibility, viz., PL is from the 
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FIG. 2: Calculated PAs from \CTX1) in pairs of (a) PPV3, and (b) PPV4 oligomers. All energies 
are relative to that of the neutral optical exciton. Arrows denote energies of absorption from the 
ground states of charged chains. 

exciton but PA is from a bound polaron-pair, is precluded by the absence of room tempera- 
ture infrared-active vibration (IRAV) in the recent PA experiment.— Unlike the CT exciton, 
the polaron-pair wavefunction is not affected by interchain charge-transfer and is asymmet- 
ric: i.e., a specific chain is positively charged while the other is negatively charged. The 
counterpart related by symmetry, with the charges on the two chains reversed, is missing in 
the polaron-pair, and hence PA even from the bound polaron-pair should be accompanied 
by strong IRAV. Although strong IRAV has been reported by one research group,— this 
results has not been reproduced by any other group. 



A key weakness of our theory in Reference 
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was that the PA spectrum could be calcu- 



lated only for a limited energy range. This is remedied in the present work. Our calculated 
PA spectra of Fig. 2 are to be compared with the experimental PA spectra of Sheng et 
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al. (see Figures 1 and 4 in Reference Il2l). Remarkable qualitative similarities between the 
calculated and experimental spectra are found, although from the theoretical spectrum it is 
difficult to distinguish between the two possible origins of the experimental Pi absorption. 
From our calculations, the Pi absorption observed experimentally can be due to either the 
low energy charge-transfer absorption, or the absorption to the \PP1*) state. We suggest 
polarization-dependent experimental measurements to distinguish between these two possi- 
bilities: the charge-transfer absorption is polarized perpendicular to the molecular planes, 
while the absorption labeled PP1 is polarized along the chain direction in the molecular 
plane. Yet another possibility is that the experimental Pi absorption is composed of two un- 
derlying absorptions. As pointed out in the previous section, in systems with realistic chain 
lengths these absorptions will be proximate in energy and can even be very close. In a recent 
experiment performed by Vardeny et al. on films of PPV derivatives, it has been found that 
the Pi absorption splits into two distinct peaks under hydrostatic pressure. 38 Whether or 
not these two PAs have different origins, as might be expected from our calculations, or 
they arise from the splitting of the same absorption under the experimental conditions, is 
currently not clear, and further experimental work is necessary to clarify this issue. 

In conclusion, we have extended the Mulliken theory of charge transfer to the case of ex- 
cited state charge transfer, to develop a theory of the CT exciton in thin films of 7r-conjugated 
polymers. The central feature of our theory is that PA in films can occur from both the 
neutral exciton component of the CT exciton as well as from the charged bound polaron- 
pair component. PAs from the exciton components occur at nearly the same wavelengths 
as those of the single chains that occur in dilute solutions. PAs from the polaron-pair com- 
ponent occur at nearly the wavelengths where ground state absorptions of charged chains 
occur. Interestingly, the absorptions from the polaron-pair component of the CT exciton are 
vanishingly weak in SCI calculations and become visible only upon including higher order 
CI. In addition to these features, the PA spectrum can exhibit a low energy charge-transfer 
absorption. Whether or not the low energy Pi absorption seen in PA measurements^ is a su- 
perposition of two fundamentally different absorptions is currently not understood, and will 
require further experimental work. There is no such ambiguity in the case of the absorption 
we have labeled PP2: this is certainly the origin of the P2 absorption in the experiment.— 
Note, however, that the relative energies of PAs labeled PP2 and PA' 2 in real systems can be 
slightly different from those predicted in Fig. 2, since electron-phonon interactions, ignored 
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in our calculations, can have a strong effect on the energy of the PP2 absorption. The 
relative energies of these experimental PAs can perhaps also be dependent on chain length, 
and hence may be more strongly conformation-dependent than the lower energy PAs. This 
may explain the difference in the conclusions regarding the energy location of P2 between 

Work is currently in progress to include electron-phonon interactions 



references 
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and 
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in our PA calculations to probe this aspect further. 
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